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O ' Abstract: NLO QCD fits using the DGLAP formalism have been made to the high precision ZEUS ■ ■ 



e+p reduced cross-section data and to fixed target structure function data, in order to determine 
parton distribution functions and the value of as(M|), taking full account of correlated systematic 

M 



errors. 



where Y± = 1 it (1 — y)^. The structure functions F2 and xF^ are directly related to non-singlet 
and singlet quark distributions, and their dependence, or scaling violation, is predicted by pQCD. 
At ^ lOOOGeV'^ F2 dominates the charged lepton-hadron cross-section and for x ^ the gluon 
contribution dominates the evolution of F2, such that ZEUS data provide crucial information on 
quark and gluon distributions. (Schematically, F2 ~ xq, dF2/dlnQ'^ ~ agPqgXg). 
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T— I ■ The new high precision ZEUS data on the neutral current e^p reduced cross section are fit to 

' the predictions of NLO QCD using the DGLAP equations, in order to determine parton distribution j^y* 
00 functions and the value of as(M^). The increased range (6.3 x 10^^ < x < 0.65, 2.7 < < 
^ SOOOOGey^) and precision (systematic errors ~ 3%, statistical errors <, 1%, for 2 < < SOOGeF^) of 
T— I . the ZEUS data allow a much improved determination of the gluon and sea distributions compared to ^ 
' our previous work [^. In recent years more emphasis has been placed on estimating errors on extracted IN3 
parton distributions. We present parton distributions including full accounting for uncertainties from ^ 
O 1' experimental correlated systematic errors. We are also able to measure the value of Qs(M^) taking full ^ — ^ 

O ^ account of the correlations between the shape of the parton distribution functions and as- Full details | ^ 

^ of the analysis are given in [^. 

^ i This analysis is performed within the conventional paradigm of leading twist, NLO QCD, with the 
• w^ . renormalisation and factorization scales chosen to be Q^. The heavy quark production scheme is the 
^ , general mass variable flavour number scheme of Roberts and Thorne |Q] 

\ The kinematics of lepton hadron scattering is described in terms of the variables Q^, the invariant mass 
of the exchanged vector boson, Bjorken x, the fraction of the momentum of the incoming nucleon taken 
by the struck quark (in the quark-parton model), and y which measures the energy transfer between 
the lepton and hadron systems. The differential cross-section for the process is given in terms of the 
structure functions by 

'^''^ [Y+F2{x,Q^)-y^Fdx,Q')-Y.xF,{x,Q^, 
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We have performed a global fit of ZEUS and fixed target DIS data. The fixed target data sets used 
are those with precision data for which full information on the correlated systematic errors is available 
(NMC, E665, BCDMS muon induced F2 data on proton and deuterium targets and CCFR u xF^ 
data on an Fe target [||). These data are used to gain information on the valence quark distributions 
and the flavour composition of the sea, and to constrain the fits at high x. However, our focus is on 
the additional information to be gained from the new precision ZEUS data, particularly on the gluon 
and quark densities at low x and on the value of as(M|). 

In the standard fit the following cuts are made on the ZEUS and the fixed target data: (i) W"^ > 
20 GeV^ to reduce the sensitivity to target mass and higher twist contributions which become important 
at high X and low Q^; (ii) > 2.5 GeV^ to remain in the kinematic region where perturbative QCD 
should be applicable. 

The QCD predictions for the structure functions needed to construct the reduced cross-section are 
obtained by solving the DGLAP evolution equations at NLO in the MS scheme. These equations 
yield the quark and gluon momentum distributions (and thus the structure functions) at all values of 
provided they are input as functions of x at some input scale Qq. The parton distribution functions 
(PDFs) for u valence, d valence, total sea, gluon and the difference between the d and u contributions 
to the sea, are each parametrized by the form 

pixP^l-xfil +P5X) 



at Qq = 7GeV^. Thus the flavour structure of the light quark sea allows for the violation of the 
Gottfried sum rule. We also impose a suppression of the strange sea of a factor of 2 at Qq, consistent 
with neutrino induced dimuon data from CCFR. 

The parameters pi — ps are constrained to impose the momentum sum-rule and the number sum-rules 
on the valence distributions. The gluon distribution has p^ = 0, because non zero values have minimal 
effect on the of the fit, and because this choice constrains the high x gluon to be positive without the 
need for penalty terms. There are 11 free parameters in the standard fit when the strong coupling 
constant is fixed to as(M^) = 0.118 and 12 free parameters when as(M^) is determined by the fit. 
Full account has been taken of correlated experimental systematic errors as follows. The definition of 
the is 

2 [-Fi(p,s) - Fi(meas))]^ ^ 2 

^ =2^ — (^2 , ^2 ) — + !L"A 

j K^stat ^ 'June) ^ 

where 

_ F,{NLOQCD)ip) 

The latter equation modifies the NLOQCD prediction for the measured quantity F (structure function 
or reduced cross-section), which is a function of the PDF parameters to include the effect of the 
correlated systematic errors. A*^** is the systematic error on data point i due to source A and sx 
are systematic error parameters. These systematic errors are assumed to represent a one standard 
deviation error. The systematic error parameters are fixed to zero for the fit, but allowed to vary for 
error analysis, such that both the error matrices 

1 ^ 1 5^x' 



2 6pj6pk 2 6pj6sx 

are evaluated, and the systematic covariance matrix is given by V^^ = M^^CC^ M^^ , so that the 
total covariance matrix is given by = VP + VP', where VP = M'^. This method is equivalent 
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to offsetting each systematic parameter by ±1, reminimizing and adding the resulting deviations, Ap, 
from the parameters determined in the standard fit in quadrature [^. Thus it produces a conservative 
estimate of the effect of systematic errors. Normalization errors are included as correlated systematic 
errors. In total 71 independent sources of systematic error are included. 

A good description of the structure function and reduced cross-section data over the whole range of 
from 2.5 to 30000 GeV^ is obtained. The quality of the fit to the new ZEUS NC reduced cross-section 
and F2 data is shown in The quality of the fit to all the data sets may be judged from the x^- 
Adding the statistical and systematic errors in quadrature gives a total P^r data point of 0.95 for 
1263 data points and 11 free parameters. 

The parton distributions extracted from the fit are shown in Fig. 
The ZEUS data are crucial in determining the gluon 
and the sea distributions so we consider these in more 
detail. In Fig | we show the extracted gluon distribu- 
tion in several bins. The uncertainty in these dis- 
tributions is further illustrated in Fig || where the ra- 
tio of the error bands to the central value is shown. In 
these figures the innermost error bands show the sta- 
tistical and uncorrelated systematic error, the middle 
error bands show the total experimental error includ- 
ing correlated systematic errors and the outer error 
bands show the additional uncertainty coming from 
variation of the strong coupling constant as(M^), 
which is taken into account with full correlations by 
allowing as{M^) to be a parameter of the fit. 
For ^ 5 GeV^, the gluon distribution rises steeply 
at low X . It is well determined to within ~ 10% for 
> 20 GeV^ 10""^ <x < 10"^ and its uncertainty 
decreases as evolves upwards. Considerable uncer- 
tainty remains for x > 0.1, where one needs information from Tevatron jet studies or prompt photon 
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Figure 1: The glue, sea, u and d valence distri- 
butions extracted from the standard ZEUS fit, at 
— lOGeV^, with the error bands resulting from 
statistical, uncorrelated and correlated systematic 
errors. 
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Figure 2: The gluon distribution from the stan- 
dard ZEUS NLOQCD fit. 



Figure 3: The uncertainties on the gluon distribu- 
tion from the standard ZEUS NLOQCD fit, shown 
as the ratio of the error bands to the central value. 
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data. At low the gluon shape is flatter at low x, and if the fit is extrapolated back to = 1 GeV^ 
the gluon shape becomes valence like, with low x values that are negative within errors. (This is also 
true of the structure function Fl which is closely related to the gluon.) The new precision ZEUS data 
collected in 96/97 shows this tendency even more strongly than the previous 94/5 data 

In Fig I we show the extracted sea distributions in several bins. The uncertainty in these distri- 
butions is is less than ~ 5% for ^2.5 GeV^ and 10~^ < a; < 10~^. It can be seen that even at the 
smallest ~ 1 GeV^ the sea distribution is still rising slightly at small x. In the same Q ,x region 
the gluon distribution is falling and is consistent with zero. Thus the rise in the sea distribution cannot 
be driven by the gluon g ^ qq splitting in this region. 

We also make a fit using only ZEUS data, including 
the charged current e~^p reduced cross-section data [^] 
and the neutral and charged current e~p data |^ , 0] . 
This high data is very well described by the stan- 
dard fit. However now we use this data to constrain 
the valence distributions. The dy valence distribution 
for this fit is shown in Fig ^ where the dy distribu- 
tion for the standard fit, including fixed target data, 
is shown for comparison. Clearly the standard fit 
gives a more precise determination, but the high 
ZEUS reduced cross-section data is able to provide 
some constraint on the valence distributions. We can 
see that ZEUS data prefer a larger d valence distri- 
bution at high X than the standard fit. 

NLO DGLAP QCD fits give good descriptions of F2 
data down to values in the range 1 - 2 GeV^ but 
such fits assume the validity of the NLO DGLAP 
QCD formalism even for low (where as is becom- 
ing larger such that NNLO corrections are increas- 
ingly important) and very low x (where ln{l/x) re- 
summation terms could be important). The fits also 
ignore high density and non-perturbative effects (see 
Refs. for a discussion of these effects). To inves- 
tigate if there is a technical low limit to the NLO 
QCD fit (in the sense that the fit fails to converge 
or gives a very bad x^)i we extrapolate our standard 
fit into the region covered by the very low ZEUS 
BPT PI data. In Fig | the main ZEUS data sample 
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Figure 4: The sea distribution from the standard 
ZEUS NLOQCD fit ^eus 
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Figure 5: The xdy distribution from the NLOQCD 
fit to ZEUS data alone. 



and the ZEUS BPT data are shown in very low 
bins compared to the fit predictions. 

We see that NLO DGLAP fit is unable to describe the very precise BPT data for < o.65 GeV^ 
even when the full error bands on the fit due to the correlated systematic errors are included. This 
remains true even if these data are included in the fit: they produce very poor fits, with per data 
point of 3.9 or more. The predictions for Fl for these low values are also significantly negative for 
^ 0.8 GeV'^ indicating that the NLO DGLAP formlism cannot be valid in this region. 



CD 

o 
o 



-4- 



International Europhysics Conference on HEP 



Amanda Cooper-Sarkar 



The value of asiM^) is measured by allowing it to be 
a parameter of the fit. We obtain, 

as{Mz) = 0.1172±0.0008(stat+nncorr)±0.0054(corr) 

where the first error represents statistical and uncor- 
related systematic errors and the second error rep- 
resents the correlated systematic errors from all the 
contributing experiments, including their normaliza- 
tion errors. The contribution of normalization errors 
is as large as half the correlated systematic error. 
The renormalization and factorization scales used in 
the fit may be varied from — > 2Q^, as a way 
of estimating the importance of NNLO and ln{l/x) 
terms. This produces variations in as{M^) of 0.004. 
Within the paradigm of NLO DGLAP the contribu- 
tion of the correlated systematic errors is the most 
significant source of uncertainty. Variation of analy- 
sis choices, such as the form of the parametrization at 
Qq, the value of Qq, the minimum of data enter- 
ing the fit, do not produce a significant model error. 
A more significant choice is that of the heavy quark 
production scheme. Repeating the standard fit using 
the FFN scheme or the ZMVFN scheme produces a 
variation in q;s(M|) of ±0.001. 

The well known correlation of the value of as(M^) to the gluon shape is accounted for by performing 
a simultaneous fit for as(M^) and the PDF parameters. The largest uncertainty in the gluon shape is 
now at large x. Tevatron jet data suggests a larger gluon at high x than that found in our standard fit, 
and if we use such data to constrain our fit we obtain as(M^) ~ 0.121, a variation within the quoted 
error. 

Thus, when a full accounting of correlated systematic errors is included, there is no contradiction 
between the lower values of as(M^) originally determined from fixed target DIS data, e.g. as(M^) = 
0.113 from BCDMS and the higher values favoured by LEP data a^(M|) ~ 0.120 [|. 



2 
1.5 

1 

0.5 


2 

1.5 

1 

0.5 


2 

1.5 

1 

0.5 


2 

1.5 
1 

0.5 




Q^=0.3 GeV^ 


Q^=0.4 GeV^ 


Q^=0.5 GeV- 


\ 

- 


- 


' *^ 

\ 

- 


Q^=0.585 GeV^ 


Q^=0.65 GeV^ 


Q^=0.8 GeV^ 


- V-., 

- , 1 ,1 ,1 ,1 ,1 






^ ^ 

- , 1 ,1 ,1 ,1 ,1 , 


- , 1 ,1 ,1 ,1 ,1 , 1 


Q^=1.5 GeV^ 


Q^=2.7 GeV^ 


Q^=3.5 GeV^ 




- , 1 ,1 ,1 ,1 ,1 




Q^=4.5 GeV^ 


C^=6.5 GeV^ 

- \ 
\ 


-5 -3 -1 
10 10 10 1 

— ZEUS NLO-QCD Fit 
(Prel. 12001 

• ZEUsyfi/? 

• ZEUS BPT 97 

• BCDMS 



10 10 10 110 



10 10 ''l 



Figure 6: F2 data down to the very low (includ- 
ing BPT data) compared to the standard NLOQCD 
fit backward extrapolated. 
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